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PROPAGATION OF A LASER BEAM IN A PLASMA 


Jon M. Chapman 
Jiraxr Kevorkian 
Loren C. Steinhauer 
Juris Vaqners 

ABSTRACT 

This paper shows that for a nonabsorbing medium with a 
prescribed index of refraction, the effects of beam stability, 
line focusing, and bear distortion can all be predicted from 
simple ray optics. When the paraxial approximation is used, 
diffraction effects are examined for Gaussian, Lorentzian, 
and square beams. Most importantly, it is shown that for a 
Gaussian beam, diffraction effects can be included simply by 
adding imaginary solutions to the paraxial ray equations. 

A .'so presented are several procedures to extend the paraxial 
approximation so that the solution will have a domain of 
validity of greater extent. 
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PROPAGATION OF A LASER BEAM IN A PLASMA 


1.) Int roductio n 

The process of heating u plasma in a solenoid using a 
laser beam is difficult to describe analytically when one 
wants to examine in a self-consistent manner, the coupled 
behuvior of plasm"' heating and beam propagation. Some insight 
into the beam behavior may be obtained however, when the ion 
density is arbitrarily prescribed. Consequently, heating has 
only been partially accounted for by allowing the ion density 
to change without an accompanyinq loss of energy from the 
laser beam. This also implies that in the plasma, only the 
electron 'fluid' is excitable. 

The impressed magnetic field along the solenoid is assumed 
to be zero since its effect upon beam propagation is small. A 
constant magnetic field in the direction of beam propagation 
effects right and left circularly polarized light differently 
so that Faraday rotation is observed for a plane polarized 
wave. However, for frequencies much larger than the cyclotron 
and plasma frequencies, the plane of polarization is essentially 
fixed and the impressed magnetic field can be ignored (1). In 
the case of interest, the ratio of the cyclotron frequency to 
the exciting laser frequency is for a magnetic field 

of 100 kG and a laser wavelength of 10. € u . 
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The cold plasma approximat ion requires that 





where is the speed of liqht, A is the Boltzmann constant, 

and " IS the mass of an electron. electron temperctures 

of 1C <i V the above ratio is 0 ( I ^ 

The above assumptions motivates the study of beam propaqa- 
tion through a cold, col 1 is lonless and field*free plasma. The 
dimensionless parameters that are pertinent to the discription 
of the interaction are obtained by nondimensiona 1 izinq Maxwell's 
equations and the momentum equation for an electron fluid. 

These parameters are then determined to be , and « ; 

the first measures a typical electron (ion) density to the 

critical electron density k , the second measures the elec- 

€(? 

tron quiver velocity to the speed fo liqht C , and the 

third parameter compares the laser wavelength • to the 
cylinder circumference .'ti . The parameter n is also given 
through the ratio of the plasma frequency to the exciting 

laser frequency u. . We thus have 


Here 
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«*U’C‘.ric field is Jetcrninel i r')r the beam intensity 

( ahu are {physical constints). 

If the laser intensity is of the order of . ‘‘ .’itij — ‘ 
jr.il the wavelength of the laser bear is » then *' • 

. * . Also, tor ion densities of ,* ** •- * and a cylinder 

rail jS of I we find that » ’ and 

The effect that V nas on beam propaq.>tion may be ander- 
srxjo uy considcrinq a plane wave traveling through a honoaeneous, 
col"f, and field-free plasma. This analysis is straight forward 
ana will not be presented here. We find that the transverse 
components of the electric and maqnet ic fields arc unchamed 
to and that sinusoidal charge separation occurs with 

a magnitude at twice tht- exciting frequency. Conse- 

quently, Its t ffects upon b«.-ar propagation are small when 
" » ; ’ and will b«.' lijnored here. 

The index of refraction »i is then easily obtained using 
the above simplifications. We find that 

' ( 1 . I ) 

where is the nondincnsionalized ion density. In this 

paper wc are going to analyze the case when k satisfies the 
relation 

fr.) - (r. ' (1.2) 
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Here k(ti and are determined fron the ion density 

and are arbitrary functions of the parameter z , the axial 
distance along the cylinder nondimensionalized by the cylinder's 
radius ^ . The coordinate r is the nondimensiona 1 radial 
distance from the cylinder axis. Prom Eq. (1.2) we see that 
the ion density is always parabolic but is allowed to vary 
with axial distance. Time variation of the ion density involv- 
ing the hydrodynamic time scale will be ignored in this 
presentation . 

The parabolic assumption is useful when the beam boundary 
IS contained well inside a favorable density profile so that 
'beam trapping* occurs. This approximation has been discussed 
by Steinhauer (2) and Mani, et al. (3] and is coupled with the 
additional approximation of extending the effective boundary 
conditions at the cylinder wall to infinity. 

The preceding discussion thus motivates the study of the 
following reduced wave equation for the transverse component 
of the polarized laser beam. 





} E 


0 


(1.3) 


The symbol represents the Laplacian in cylindrical coor- 

dinates. Time has been nondimensionalized w.r.t. the exciting 
frequency Ji so that is understood to include the factor 
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xs obtained when 


An asymptotic representation for £ 
the intensity of the incident beam is neqliqible in reqions 
wnerc the ion density is close to or qreater than the critical 
ion density (this condition ir. fact motivates the parabolic 
approximat ion . This asymptotic representation is easily 
obtained from the exact sciut’jn of Eq. (1.3) when and .. 
in Eq . (1.2) are constants. This procedure will be developed 
in Section 2. Alternatively, the asymptotic representation 
for can be obtained directly from the differential Eq. (1.3) 

ari will be presented in Section 3. 

When k and u are constants the usual asymptotic repre- 
sentation predicts a periodic structure for the beam propaqa- 
tion (as does the standard paraxial approximation) and thus 
does not include 'washout* effects of the beam's intensity 
distribution. The term washout means that the oriqin.il inten- 
sity distribution (at r. * C ) can never be reqenorated for 
any value of z and that it continually chanqes for increis- 
ingly larger values of the variable r. . Ray optics predicts 
this distortion effect from the result that rays with different 
initial displacements from the axis have different focal lengths 
in the parabolic focusing medium [4]. It is because of this 
property that line focuses appear with monotonical ly increas- 
ing lengths. 
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In Section 2 wi> will exairine the 'beaa boundary* describ- 
ing the propagation of Gaussian beams. Susw results are also 
given for beams that have nonGaussian distributions; in parti- 
cular we will look at a beam that is initially Lorentzian and 
a beam that is initially S'|Uure. 

In Section 3, of prime importance, is the development of 
a less complicated equation which describes the beam boundary 
t.'a^ectories for Gaussian beams. Analysis here will be mostly 
Jevoted to the case when and ^ ire *sl^ly* varying func- 
tions of z . However, other special situations will be 
examined. In particular, we will see that instabilities can 
occur even in a 'favorable’ density profile. 

In Section 4 we will see how the normally neglu~ted wash- 
out effects of the bean’s intensity distribution may be 
included in the asymptotic representation. 

Section S includes for completeness, the analysiM of 
the ray optics equations when • and .. are slowly varying 
functions. The idea here is that a uniformly valid solution 
can be obtained in reqions away from the line focuses. Con- 
sequently, washout effects have been obtained at the expense 
of neglecting diffraction. 

Section 6 considers higher order terms for a Gaussian 
beam. The small secular effects of beam distortion and line 
focusing are then examined. 
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2 . ) The Aayinptot ic Development from the Exact Solution; Fixed 
Parabolic Profiles 

The solution for the transverse component of the electric 
field can be represented through an infinite sum of 

Laquerre polynomials. To accomplish this introduce t*ie complex 
field amplitude expressed in transformed variables o and C 
sc that 

tKZ 

« c \l)(p,Q,rJ 

where 

0 • ur* 





( 2 . 1 ) 


here > and u are constants. 

The differential equation for ii> , is obtained from 

r.q . (1.3). Th us 
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* 'K 
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kLi . 


f.f. 


( 2 . 2 ) 


.Separation of variables and the physical requirement of bounded- 
ness of ’ll as () approaches infinity gives rije to the expo- 
nential Laquerre polynomials as eigenfunctions of Eq. (2.2). 

In terms of these eigenfunctions i|/ can be expre.^sed as 

y r fr,t) l (^) (2.3) 

o " " 
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ki 


where Laguerre polynomial of order n and 

T^a.t) ~ I L^(p) dp (2.4) 

*0 

Applying the transform defined by Eq. (2.4) to Eq. (2.2) gives 
the differential equation for • Thus 

^ T •• -h 2i r • - 4(n4-h r ~ 0 (2.5) 

k* n n 2 k 

The prime denotes differentiation w.r.t. ^ . 

The roots c. of the characteristic equation corres- 

n \,2 

ponding to Eq. (2.5) are 


n \,2 


- ^ r; ♦ /7 

UE* 






( 2 . 6 ) 


The choice of the minus or plus sign corresponds to right or 
left propagating waves respectively. Since we are interested 
in beam propagation down a solenoid the second root can be 
ignored. 

Several observations considering the exact solution may 
now be made: 

(1) Different eigenfunction solutions have different 
periods in C so that a solution representing a sum 
of such terms cannot be strictly periodic. 
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(2) If the initial distribution c«r be represented by a 
sum of eigenfunctions of relatively low order then the 
solution will be nearly periodic. This will be the 
basis of the asymptotic analysis. 

(3) To each eigenfunction there is an associated bundle 
of rays located at a particular initial radius from the 
axis. This is due to the fact that ray optics predicts 
different periods for rays with different initial dis- 
placements. For larger values of n we obtain larger 
initial displacements. 

(4) By the analogy of (3) we see that a damped high order 
eigenfunction solution corresponds to a ray entering the 
medium with an initial off axis displacement large enough 
so that the ion density exceeds the critical density. 

We now s^.ek an asymptotic expression for tl> , r , i in 

the following way. Let the infinite sum in Eq. (2.3) be 

replaced by a finite sum of the first A/ terms. If is 

large enough this sum will converge to the exact solution. 

Now let c be small enough so that 


4 (N ■¥ < 1 . 


The square root in Eq . (2.6) can now be expanded. We find that 


u » i(?.n+l) 0(r^ ) V n < N 

n 1 


(2.7) 


The expansion for is then represented by the series 


r fr„c) = ^ * 

n K 


} ( 2 . 8 ) 


We also obtain a corresponding series for ; . Thus 
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i,c) 


-C/Z I ^-ilZnHH 

*' n FI 

0 

» **•*’ 

0 


The coefficients c are to be determined from the initial 

FI 

conditions. 

Since S is a measure of the degree of accuracy by which 
one can replace an infinite sum by a finite sum, we expect 
that the first sum in Eq. (2.9) may be considered in the 
limit as S ■* <*> . 

Let be defined by this infinite sum. Thus 




-p/2 




LniPJ 


( 2 . 10 ) 


Consequently, the solution for can be expressed 

asymptotically by the series 


>M'P, 


i- 




♦ . . . 


( 2 . 11 ) 


We sec that the representation is useful when both the secular 
term ye*? and the curvature term | | arr small. When 

the latter condition is fulfilled a 'multi-v&riable* procedure 
[51 can be used to extend the domain of validity of the asymp- 
totic representation. 
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It is evident that in periodic in and that 

downstream 'washout-' effects are not included. If uc^L <<■ 1 
where L is ♦he extent of the domain of ^ that is of inter- 
est, then the distortion effects can be considered by includ- 
inq the next term in Lq . (2.11). However, if ‘ i » O.’.' 
then one must resert to the multi-variable procedure of 
Section 4. 

To obtain in a 'summed' form we must consider the 

coef : ’ cients given by the relation 




r 

j 


-p/i 


/, U. 


vfr. 


( 2 . 12 ) 


If we introduce this expression into Eq. (2.10) and inter- 
change summation and integration, the following expression 
can be obtained 




I 


£ u/.*; cot ' r’ 
■ ’’ sin <•'. j 



< nJ 


I sin ■ I 


(2.13) 


We w*ll consider three cases of beam propagat ion th.it involve 
an infinite sum of Laguerre polynomials. These cases may bo 
designated by their initial field distribution. The designa- 
tions arc 

■ 

i) Gaussian 
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ii) Lorentzian 




iii) Square 




-v/2p 

« 

1 

( I 0 < ►.-/* 

0 P > Uki* 


(2 


Here , s , and d are constants defined through the 
given intensity distribution and v is consid* ed to be 
sma 1 1 . 


When the beam is initially Gaussian the integral in 
Eq. (2.13) can be easily evaluated. We find that 





h(V 


cos C ♦ i 





(2 


We may also represent the complex function h(t.) by its 
amplitude a(^)/a^ and phase so that 


h(U 



and 


a • a(0) 
o 


(2 


In terms cf these real functions we find that 





exp 


ftp di 

^2a 



12 


.14} 


.15) 


.16) 


k 


with 


and 



(2.i7) 


Wo see from Eqs. (2.17) that i ' ) is the so-called 'beam 
boundary* because it qives the value of .•* at which the 
intensity is do%m by a factor of . as compared to its 
value alonq the axis. Moreover, ve find that the ray paths 
for are simply described by some constant times . 

These results concerninq Gaussian beams have been obtained 
previously by various authors [3,6]. 

The inteqral in Eq. (2.13) is harder to evaluate for 
non-Gaussian beams. There may be special circumstances 
however, when one miqht use to advantaqe the complex con- 
volution theorem concerninq J.^placf transforms. To qiven 
an example of this procedure we will consider the Lorentizian 
beam. 

Let /V.?; » L|/| denote the Laplace transform of a 
function f(^) . For the Lorentzian beam we have 



(2.18) 
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H«r« E (*J is ths sxponsntl:*! Intsgral written as 




(2.19) 


We are interested in the case when 


8 » 2 (\f-i cot 


( 2 . 20 ) 


Since 





fyH» 

I 

Y-t« 


F (a~{i)) du 
2 


( 2 . 21 ) 


the integral in Eq. (2.13) can be expressed in convoluted form. 
Let 1(a) represent this integral. We find that the 
contour of 1(a) may be closed to the right so that 

b 

- FT I ^ [i '2- “> 


There is an essential singularity at the point u » 5 . 
Consequently, the integral can be evaluated by Cauchy's 
integral theorem of complex analysis to give 

- f ?sm fr [s <2-“' 

Needless to say, this form is also difficult to interpret. 
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It may be of use however, in that it may converge faster 
than the corresponding infinite s€‘ries of Laguerre polyno- 
mials. Of course, Eq. (2.22) is of interest in itself in 
that it represents an alternate expression for the integral 
in Eq. (2.13) and that powerful techniques of analysis are 
available to examine its behavior. 

We next as)c if Lhere are any special values for C in 
which the integral of Eq. (2.13) can be simply evaluated. 

We see that the focal point f, = u/3 is such a point and we 
obtain the following results: 

i) Lorentzian *^ 0 ^*^* ^ 

ii) Square \ii^(p,^) * - id /u/p J ^ (d/up) (2.24) 

Here K is the modified Bessel function of order zero and 
0 

fJ is the standard Bessel function of order one. Diffrac- 

t 

tion rings at the focal point are apparent for the square 

Iseam because of the zeros of J . However, when (. * it 

1 

these fade away and the solution regenerates its initial 
distribution. This is just the periodic behavior of 

An alternate approach for the description of beam propa- 
gation can be accomplished with use of higher order Gaussian 
modes since any beam can be considered to be an infinite 
sum of such modes. The advantage of such an approach is 
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that the ray paths of each node can be easily detemined; 
the disadvantage is that the solution still involves an 
infinite sun. 
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3 . ) Changing Parabolic Profile \ 

The asymptotic representation fot can be obtained 

directly from the differential equation (1.3) when k and u 

in Eq. (1.2) are arbitrary functions of the scaled variable 

n « Ls . Thus 
1 

k'r.,- '■ * 

and n< r.,( ) * utr. (3.1) 

I 

Physically, this has the interpretation that the ion density 
is allowed to vary on an axial scale length comparable to 
or slower than the focal length of the parabolic focusing 
med i urn . 

We now obtain the equation for the complex field ampli- 
tude ^ through the transformations 

(r.) 

k — ■■ 

i * v(r,:: ,a < 

I 

where » l.': 

i 

•ir. 

and - 7 - * k(:: } (3.2) 

'is 1 

The differential equation for is then given as 
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We My no» n« 9 l*ct th« last tara in Bq. (3.3) %#hanavar the 

a 2 ^ 

curvatura tjra If— rl is smII. Froai our pravious analysis 
we sea that this is just the paraxial approxination and con- 
sequently washout and line focusing affects are neglected. 

When the beam is initially Gaussian, then the equation 
for the beam boundary a(»^) as expressed by Tien, at al. 

(6) can be obtained through the following transformation. 
Assume 


♦ Tr.a.c; 



(3.4) 


Substitution of this equation into Eq. (3.3) gives the 
differential equations for the transformed functions a(u^) 
and ) Thus 

I 


k^iz 




(3.5) 


k(z ) 4^ • 

I 


2* ^ 


(3.6) 
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These equations can be written in simpler form if we 
introduce the variable a * where 


~3m * k(z ) 

I I 


a (0) • 0 

I 


(3.7) 

(3.8) 


!• 



With this transformation Eqs. (3.5) and (3.6) become 


d'a 

dz~^ 


+ )a - • 0 

I a* 


(3.9) 


1 


(3.10) 


The function u(fi ) is to be considered as an implicit 

1 

function of z * obtained through use of Eq. (3.7). 

Obviously the analysis of Eq. (3.5) or (3.9) is diffi> 
cult because of the nonlinear term. One can obtain however, 
a linear equation that describes the beam boundary trajec- 
tories. To do this we assume that 




€) 


./k(W 1 

“ /k(z ) h(z *) 

I i 


tr* (1 (ih \ 

2 dz~f^> 
e I 


(3.11) 


The differential equation for h(z^*) is then easily 
determined to be 


(3.12) 

\ 

This is exactly Eq. (3.9) without the nonlinear term. 

Comparing Eqs. (3.11) and (3.4) gives the following 
relationships 


h(z *) 

I 


a(z *) 


oToT 



(3.13) 
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(3.14) 


1 dh 1 dx ^ i 

JT ar? • j asT? ♦ ? 

I > 

a( M *) 

From Eq. (3.13) we see that a(0) is just the amplitude 
of the complex function k(M^*) and is the phase. 

To determine the solution of Eq. (3.12), initial condi- 
tions must be specified. These can be obtained from Eqs. 
(3.13) and (3.14) when the initial conditions on the beam 
boundary are )cno%m. When a(0) » <Jq * (Oj “ and 

I 

^(0) " 0 ve obtain 


h(0) - 1 


when 


dh 


( 0 ) 



k • k(0) 
o 


(3.15) 


If the initial conditions on k were arbitrarily speci- 
fied to be real then Eq. (3.12) describes the trajectories of 
the beam's center (6]. Consequently, by allowing complex 
values for h(z^*) , we are able to describe the trajectories 
of the beam's boundary. Moreover, since the real and imaginary 
parts of h(s^'*') satisfy the same equation, stability or 
instability of the beam's center occurs together with t..e 
steUsility or instability of the beam's boundary. 
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The solution of ttq. (3.12) is readily determined when 
k and u are slowly varying. We measure this slowness through 


the small parameter v so that 


where 


•:(r. ,'j) 

i 




3 • ^r. 

. I 


( 3. 16) 


We now represent h.s *,\f asymptot ical ly and assume this 

to iepemi open two axial scales Z and so that 


h(: 


.v> 


^ ' I" •• 1 
j 


♦ jn 




iz.z ; 

7 


(3.17) 


where 7. is defined by 
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The derivativ** term in F.q. (3.12) can bo obtained through 
the relation 


1 


lz~* 


I 



r. 

7 



( 3. 19) 


Notice must also bo made of initial conditions and we obtain 
these from (3.15) . 
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standard multiple acale techniques then yield 


2 0 0 0 0 


(3.20) 


Observe that the real part of this solutior is just %irhat 
one would obtain from ray optics for a paraxial ray. The 
addition of an imaginary term satisfying special initial 
conditions then improves the optics approximat ion by includ- 
ing the effects of diffraction. We h^ve also remarked that 
a uniformly valid solution could be obtained by extending 
the paraxial approximatioi.. If thus seems feasible that this 
may be accomplished by allowing complex rays where the real 
part satisfies the exact ray optics equation and where the 
imaginary part satisfies special initial conditions. This 
has not yet been verified. 

A case exhibiting various regions of stability and 
instability can be examined by studying Mathieu's equation. 




*■ 


A * 2W cos ZX)h 


0 


(3.21) 


Eq. (3.12) can be placed in this form when k is constant a 
and varies sinusoidally. When A is zero the medium 

exhibits periodic focusing and the solution is first stable 
for values of « less than 0.8€ (6). Regions of stability 
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and instability alsc occur for positive A and are again 
determined by the values of Q . Consequently, instabilities 
can occur even in a 'favorable' density profile. 

The case when k has small sinusoidal variations and 
u is constant has also been examined by Tien, et al. (6]. 

A more general stability analysis is needed for the case 
when both k and u are varying sinusoidally w.r.t. the fast 
variable but are also allowed to vary slowly w.r.t. the 

parameter a « va 

2 I 
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4 . ) Extended Paraxial Approximation 


In Section 2 we obtained m ap'^roximation that was valid 
as long as the secular term uc’^<< 1 . If we define 
V>t through the relation 

nA'n-e^M* (4.1) 

so that An represents the normalized change in the ion 

density at the cylinder's edge compared to the axis when 

the distribution has the assumed parabolic profile, then 

the secular criterion can be expressed as nA'n cs << 1 . 

Since z has been nondimens ionali zed by a dimension on the 

order of 1 an and e ■ 0(10~*) we find that the paraxial 

1 0 kfn 

approximation is valid for distances less than — — . If 
/Cn is small, that is if the density distribution is fairly 
flat, then the approximation will be valid over the domain 
of interest for fusion reactors; whereas if An « 0(10*) 
then clearly we need to extend the domain of validity of 
the approximation. Note that if we require the second term 
in Eq. (2.11) tc- be less than 10% of the first term, then 
An » 0(10) would require us to examine che extended asymp- 
totic representation. 

The mevhod of analysis will be to introduce a new 'slow 
variable' T = ^ and to obtain the necessary functional 

dependence of the approximation on this variable by considering 
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both the Laquerre series solution and the differential equa- 
tions obtained directly from the reduced wave equation (1.3). 
For simplxr.’ity we will consider here only the case when 
and in Eq. (1.2) are constants. 

From the expansion of (2.6) we can easily obtain the 
dependence of upon t . Thus 


-p/£ r 

c ) a c 
0 



(4.2) 


When the beam is initially Guassian, it is hoped that we can 
sum the above series. This summation may be best accomplished 
however, if we examine the partial differential equations 
directly. 

To do this we introduce the t rans format ions 


i_ ♦ i ♦ «» 

* %■, Tf 


i/Cr.r.c; - ^ (4.3) 


The Eq. (2.2) is then transformed into the pair of equations 


and 



7 * Uu, - » 


(4.4) 


({ i- . i -i’. I 


l|) * 0 

o 


(4.5) 
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Let the function satisfy Eq. (4.4). When 

the beam is initially Gaussian, we can specify 
through Eq. (2.15). Separation of variables of Eq. (4.5) 
leads us to suspect that 


when 


C.T) 


f* \i>^(p,r)»0) dr^ 


(4.6) 




1 P -2: 

2 ^ J ^ 


dX 


(4.7) 


The integral (4.8) is quic)ily evaluated to be 

i _ ijl 

e ^ ^ 

- 


(4.8) 


so that Eq. (4.7) ta)ce8 the form 




-t ^/1 

<g 






'i>Q(P$^i3 0) dr'. 


(4.9) 


Powerful techniques of complex analysis may now be used 
to study the behavior of once it has been veri- 

fied that Eq. (4.9) is correct. We find that this is indeed 
so by substituting the Laguerre series representation for 
^^(P»^»0) from Eq. (2.10) into Eq. (4.9) to recover Eq. (4.2). 
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5. ) Ray Optics and Slow Variations 

A uniformly valid solution can be obtained in regions 
away from line focuses he neglecting diffraction. With 
this approximation, ray trajectories are simply obtained 
from the standard ray optics equation. Here we again con- 
sider a medium with an axially symmetric index of refraction 
and which is also alla»wed to vary slowly. We thus specify 
the index of refraction to be given by 

* k^(s) - u^(z)r^ 


5 » Oa (5.1) 

so that the slowness is defined by the small parameter v . 

In terms of previous notation we have the following identities. 

V » cv 

u * cu 


S a (5.2) 

I 

The ray equation for a medium with index of refraction 
(5.1) ir simply 


la 



’Jn 


in 


(5.3) 
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Introduce a new function 


pr*,v; 


through the relation 


P(z,m) 




(5.4) 


Then tne above ray equation becomes a system of two equations 
for r and P . 


. j 

cm' P 


(5.5) 


dz 



(5.6) 


By differentiating Eq. (5.5) w.r.t. z and using Eq. (5.1) 
tite above equations can be expressed in the following form 


d^r 

d^ 






0 


(5.7) 



Z\)[kk' - uC'r*]P*^ 


(5.8) 


where the new variable z'*’ has been introduced by 

dz*" p-in 
dz ^ 

z*’(0) • 0 (5.9) 

The prime denotes differentiation w.r.t. the parameter i . 
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Again z is interpreted to be an implicit function of 
obtained through use of Eq. (5.9). 

We now assume r(z*,v) and P(z*t\>) to have asymptotic 
expansions dependent upon two axial scales Z and z so that 

r(z*t\t) * r^(Z,z) + \>r^(Z,z) •#•... 

P(z*',\)) » P^(Z,z) ^ vP^rz,5; ^ . . . (5.10) 

when the new vaviable Z is defined by 

U(z) Z(0) * 0 (5.11) 

The derivative terms are obtained through the relation 

- DfiJ ^ ^ (5.12) 

Usual multiple scale techniques then yield 

r (Z,z) = J ^717 cos(Z-f * ) 

0 m V ^(z) 0 

PJZ,z) = k^(z) - u(z) (5.13) 

The constants r and 4)^ are determined from initial 

tn o 

conditions . 

It can be seen from Eqs. (5.13) that slow on axis 
density variations effects only the period of each ray 
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A 


and does not influence its slowly changing amplitude, 
result is in agreement with Steinhauer [ 4 ]. 


This 
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6. ) fliqher Order Considerations 


In Section 4, the multiple scale technique was used to 
extend the domain of the asymptotic representation so that 
• on) is allowable. On the other hand, if we are 
interested in the higher order terms with *■' ) then 

a particularly effective method of including the secular terms 
in the asymptotic formula can be obtained by 'straining' the 
coordinates o and C . The small effects of beam distortion 
and line focusing can then be examined under the stated res- 
triction << 1 

The asymptotic expansion of was found to be 


l/fr,’..t) - it). 


I yL 

r k 






( 2 . 11 ) 


The secular term can be removed from this expansion by use 
of the strained coordinates H and : defined through the 

reciprocal transformations 

f, * * i/p th,T.) * . . . 

I 

( n ^ f i' *1 . ( fi . 1 ) 

I 

From these equations one can determine ;; and H as functions 
of and p . Thus 

.T = f. - e *.• (p , . 

I 

h ^ P ~ f ^ t p,f.) * .. . (6.2) 
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The functions e and t are to be so chosen that the secular 

> I 

term in (2.11) is eliminated. This can be accomplished by 
substituting (6.1) into (2.11) expanding and collecting terms 
of Ofc‘) . This term must be identically zero. Thus 


e 

> 


dw 

0 . 




( 6 . 3 ) 


Pron this we see that we are at will to specify one of the 
straining coefficients. We thus set equal to zero for 
this choice allows readily interpretable results. We then 
must choose 


♦ 

** I 



(6.4) 


Performing the required differentiations on 
the beam is Gaussian gives the result 



iph* . U*L .... 

iy * Tir I— ST— 


T' 





when 


(6.5) 


We may express h'/h in terms of its real and imaginary 
p>arts through the relation 


Also 


c-' . , J , 


\m* 


a 

— e 
a 

o 


(6.6) 
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Usinq these expressions, the amplitude of can be found. 

Thus 


±_ 


• •'o' 


{exp[- ~ coaS^-f/- alni'tl 


- (t5 


(6.7) 


The beam has its maximum and minimum amplitude when >;i « >/tt 

T 

and ; * 7 - Lk*^1) respectively. At these points the intensity 

distribution determined f»"om (6.7) is identical to that which 
ignores the secular term. Consequently, the beam is undis- 
torted at these points. Line focusinn effects can also be 
^■tcrmined from (6.7). 


I 


3J 


k 
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